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We study the dynamics of neutrally buoyant particles with diameters varying in the 
range [1,45] in Kolmogorov scale units {ij) and Reynolds numbers based on Taylor scale 
{Re\) between 580 and 1050. One component of the particles' velocity is measured using 
extended Laser Dopplcr Vclocimetry at the center of a Von-Karman (VK) flow, accelera- 
tion is derived by differentiation. We find that although the particle acceleration variance 
decreases with increasing their diameter with scaling close to {D/rj)'^/^ , in agreement 
with previous observations, the characteristic time of acceleration autocorrelation in- 
creases much strongly than previously reported, and linearly with D/rj. A new analysis 
of the probability density functions of the acceleration shows smaller wings for larger 
particles; the flatness indeed decreases as also expected from the behavior of Eulerian 
pressure increments in the VK flow. We contrast our measurements with former ob- 
servations in wind-tunnel turbulent flows and numerical simulations, and discuss if the 
observed differences arise from inherent properties of the VK flow. 



1. Introduction 

Research in dynamics and transport phenomena in turbulence have recentl y bene - 
fitted from experimental tracking of f l ow tr a cers, see for instance ( Ott fc MannI (2000 ): 

La Porta efoZ] (|200l[ ): lMordant et all (|200lh : lArneodo et a^l (|2008hmschi fc Bodenschatz 

Ideally these tracers should have a size much smaller than the Kolmogorov 
length [r]) where the velocity gradients are smooth and hence their motion follows fluid 
streamlines, but experimental constrai nts have of te n lead to the use o f larger particles 
- with some bias as discusses e.g. in ( Meil ( 19961 ): Brown et al. (|2009l) ). On the other 
hand, the question of the dynamics of objects with a finite size freely advected by tur- 
bulent motions is a question on its own. Theories developed in the small particles limit 
at vanishing particle Reynolds numbers Rcp yield the celebrated Maxey-R iley-Gatignol 
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equati o n, with little count e r par t for high-j^ep situatiori s (see, however, 

( 1988 ): Lovalenti fc Bradv ( 1993 ): Loth fc DorganI ( 2009() ). A recent systematic analy- 
sis has been made in a wind tunnel {Rex = 160) using helium inflated soap bubbles 
(|Qureshi et all (|2007t )). Other studies were performed with isoden se polystyrene par - 
ticles in wat e r in a turbulent Von-Karman flow {Rex € [400. 815]) (jVoth et ali (|2002[ ): 



Brown et al. ( 20091 )). They have obtained several noteworthy results: (i) the variance 
of acceleration decreases as D~^^^, after the manner with which pressure increments 
vary with size according to the Kolmogorov scaling predictio n. The influence of v arying 
Reynolds numbers {Rex G [400, 815]) has also been studied in lBrown et all ( 2009f l show- 
ing that the variance of acceleration actually scales according to €^^^i^~^^^ {D /r])~^^^ . 
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(ii) The probability distributions functions (PDFs) of accele ration compon e nts do not 
depend on part icle sizes in the exp lo red range D/p £ [12 , 25] lOureshi et al. ( 2007 ) and 
D/i] e [0.4, 27] (jSrown et al\ (|2009l )V iBrown et aU (|2009l) note that PDFs for large par- 
ticles sizes may be slightly below the fluid particle PDF however due to systematic 
uncertainties they can not draw any fi rm conclusion. Our numerical study based on the 
Faxen Model ( Calzavarini et Hh (|2009D ) suggests instead that the PDF shape should be- 
come narrower at increasing the par ticle size for Re\ = 75. 180. This finding is questioned 
by a more recent numerical study i iHomann fc Bed ( 20091) ) based on a direct simulation 
approach via penalty methods, which finds a collapse of the PDFs for D/ij E [2, 14] at 
Rex = 32. We go beyond the above mentioned observations usin g new expe r hnenta l data 
obtained by extended Laser Doppler Velocimetry in a VK fiow (|Volk et al\ (|20086h V We 
find that: (i) There may be corrections to the (D/rj)^^^'^ scaling of acceleration vari- 
ance which comes from intermittency. (ii) The acceleration PDFs normalized by their 
variance are not independent of the particle size. By an indirect statistical analysis we 
observe that the wings of the distributions become less extended at increasing D/rj. (iii) 
The response time of the particle, as computed from the ac celeration autocor r elatio n 
function, increases much strongly than previously reported in Calzavarini et all (|2009l ). 
and linearly with D/rj. (iv) Still there are now consolidate differences between results in 
von Karman (VK) and wind-tunnel (WT) turbulence, such as different functional forms 
for the acceleration PDFs, or the trend in the acceleration correlation time with particle 
size. We argue that such discrepancies may originate from the large scale structure and 
anisotropy characterizing the VK flow. 



2. Experimental setup 

The flow is of the von Karman type, as described in lVolk et al\ (|2008al ). Water flUs a 
cylindrical container of internal diameter 15 cm, length 20 cm. It is driven by two disks of 
diameter 10 cm, fitted with blades in order to increase steering. The rotation rate is fixed 
at values up to 10 Hz. For the measurements reported here, the Taylor based Reynolds 
number reaches up to 1050 at a maximum dissipation rate e equal to 22 W/kg. The flow 
temperature is regulated to be constant and equal to 15°C at all rotation rates. The 
tracked particles have a density of 1.06, with diameters D = 30, 150,250,430,750 /xm. 
When further changing the fiow stirring, this will corresponds to D/jy € [1, 45]. 



Parti cles are tracked using the extended Laser Doppler Velocimetry introduced in lVolk et al. 



(1200830) ■ We use wide Laser beams to illuminate particles on a significant fraction of 
their path. When a particle crosses the fringes, the scattered light is modulated at a 
frequency directly proportional to the component of the velocity perpendicular to the 
fringes. As the beams arc not collimated, the inter-fringe remains constant across the 
measurement volume whose size is about 5 x 5 x lOmm^. In practice, we use a 2W 
continuous Argon laser of wavelength 514nm to impose a Almicrons inter-fringe, and 
image measurement volume on a low noise Hamamatsu photomultiplier in the case of 
the smallest (fluorescent) particles, while for larger particles, the detection is made using a 
PDA-36A photodiode from Thorlab. The output is recorded using a National Instrument 
PXI-NI5621 16bit-digitizer at rate 1 MHz. The velocity is c omputed from the ligh t scat- 
tering signal using the demodulation algorithm described in lMordant et al. (2002), with 
a time resolution of about 30 fis. The acceleration is then computed by differentiation of 
the velocity output. Because of measurement noise the signal h as to be filtered using a 
gaussian smoothing kernel with window width w as proposed in lMordant et ali (|2004ia ). 
Moments of the statistics of fluctuations of acceleration are computed for varying values 
of w and interpolated to zero filter width. 
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Table 1. Parameters of the flow. Q.: rotation rate of the disks, e dissipation rate, from the 
power consumption of the motors (accuracy of about 20%). The Taylor-based Reynolds number 
is computed as R\ = ^ l^uf^^ / tu , and ao is derived from the Heisenberg-Yaglom relation 
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3. Results 



3.1. Particle velocities 

One expect that Eulerian and Lagrangian velocity statistics coincide under ergodicity ap- 
proximation, so that the tracer particles velocities are expected to have Gaussian statis- 
tics. Our observation is that the velocity distribution is markedly sub-Gaussian as seen in 
figure [Ija) - flatness values are (2.56,2.58,2.62,2.46) for the 4 Reynolds number values 
explored in this work. Sub-Gaussian statistics for the velocity have been observed in many 
experimental set-ups, however usually less pronounced than in our case. Flatness values 
for velocities in WT flow arc closer to 3 (Bourgoin private com munication) . Here t he vo n 
Karman flow ha s a larg e scale inhomogene ity and anistropy (c.f. lMarie fc Daviaud ( 2004 ) 



Ouellette et ai\ (120061) : IVolk et al\ (l2006l) ) which may enhance the sub-gaussianity. In 



such a confined geometry indeed, the V K flow is known to have several possible configu- 
ration s of its large scale velocity profile ( Monchaux et al\ ( 2006 ): de la Torre fc Burguetel 



([20071)); each configuration may lead to Gaussian velocity fluctuations about a locally 
different mean value with an overall effect leading to a sub-Gaussian histogram. However, 
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Figure 2. One-component normalized acceleration variance {a^) / (a^) versus particle size. In 
order to be able to compare flows at varying Reynolds number Re\, the particle acceleration 
variance is normalized by the one measured with the smallest particles (tracers (T)) for which 
/3/?7 ^ 2 at all Re\. Solid line: Kolmogorov scaling {a%) oc (D/r;)^^''''. Dashed line: refinement 
including intermittency corrections (section 3.4). 




we have not observed any change in the velocity statistics when the Reynolds number is 
increased or when the the size of the particle is changed by over an order of magnitude in 
D/t]. In fact, for this fully turbulent regime, the velocity variance is equal to 30% of the 
impeller tip speed -HJa). This is a characteristic of the yon Karman driving impellers, 
and not a characteristic of the inertial particle size, c.f. iRavelet et al. ( 20081 ). We note 



that this observation is in agreement with a p rediction following Fax en argument at the 
leading order, wV^fluid = 1 - (1/100) (LiVA^) (jHomann fc Bed (|2009[ )) - where A is Tay- 



lor's microscale, giving a correction smaller than 1% at the Taylor-Reynolds numbers 
considered here. 

3.2. Particle acceleration variance 

With one component of velocity probed by the eLDV system and in a situation where the 
direction of motion is not prescribed, the first moment of the distribution of acceleration is 
null. One expects that the second moment (acceleration variance) reduces with increasing 
particle size, because the pressure forces which mainly cause the motion are averaged 
over a growing area. As shown in figure [H this is indeed observed. The evol ution of the 
accele r ation variance measur e d here is in a g reeme nt with previous studies bv lVoth et al 



( 2002| ): Qureshi et al. ( 2007 ): Brown et all (|2009l ): when normalized by the acceleration 



variance of the smallest particles (fiuid tracers, noted (a|,)), the quantity (a|,)/(a|,) 
exhibits a decrease consistent with the power law {D/ri)~^/^ for all the Reynolds numbers 
and inertial range particle sizes. Let us recall that this power-law behavior is obtained 
when one assumes that the particle acceleration scales like pressure increments over a 
length equal to the particle's diameter. In the inertial range of scales, this argument 
yields the scaling (al) oc {{SoP/Df) - D^/s-s = jr)-2/3 _ ^j^g^.^ (j^p)" ^ S'^{D) 

is the pressure n^'^-order structure function. We show later in text that one may also 
include intermittency corrections to obtain the dashed line shown in figure [5] which yields 
an improved agreement with measured data. 

3.3. Particle response time 

A characteristic time for the evolution of a particle response to evolving flow conditions 
is obtained from the acceleration auto-correlation functions. Their shape and evolution 
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Figure 3. (a) Acceleration auto-correlation functions fo r Rex = 815; ( b) evo lution of particle 
response times. The wind tunnel data is extracted from iQureshi et al\ (|2008l ): the Faxen data 
are from [Calzavarini et al\ (|2009l V For the eLDV data the symbols corresponds to increasing 
Reynolds numbers: (V) Rex = 590; (<) Rex = 815; (t>) Rex = 850; (A) Rex = 1050. 



with particle size is shown in figure [3l^a), for Rex = 815. In agreement with previous 
observations for tracers, the auto-correlation for small particles vanishes in times of the 
order of a few Kolmogorov time r,, = \/v/e. As could be expected, the response time Tp, 
as defined as the integral over time of the positive part of Caaij), increases with size, 
at any given Reynolds number. Our observation is that for a given Reynolds number, Tp 
increases linearly w\t}\ the particle diameter D for sizes larger than about IOt]. In addition, 
as shown in figure 3(b)(rcd/trianglc symbols), measurements performed at various Rex 
all line-up on the same curve, confirming that the evolution is indeed given by the relevant 
dimensionless variables, Tp/r,, = /(D/yy), i.e. when the response time is counted in units 
of the Kolmogorov time = -^/z^/e and the particle size counted in units of dissipative 
scale Tj = (z/^/e)^/^. 

The observed behavior is quite different from the prediction of point-particle (PP) 
models, for which the Stokes drag term becomes rapidly negligibl e when the particl e 
size increases, so that the response time remains that of fluid tracers Volk et al. ( 2008qI ). 
A first refinement of the PP-model is to account for size effects by averaging the flow 
fields over the area of the particle (for the estimation of drag) and over its volume (for 
ad ded mass effects); this is the essence of the Faxen-corrected model (FC) introduced 
bv I Calzavarini et al. ( 2009() . Using it, the authors have observed a variation of the par- 
ticle response time with size: it increases by almost a factor of 2 when the size of the 
particle increases from = 2?7 to up to I? = S2r]. Th is flnding was i n gene ral agreement 
with experimental measurement in a wind tunnel by iQureshi et al\ (|2008f ) . As shown in 
figure 3(a), our observations in a von Karman flow show a much steeper increase: the re- 
sponse time of the particles is about four times that of the tracers when the diameter has 
grown to 3277, and the variation is roughly linear. However, our data are not sufficiently 
extended to exclude the scaling Tp/r,, = {D/rff'/'^, which comes from assuming that the 
response time varies as the eddy turnover time at the scale of the particle. 



3.4. Particle acceleration probability density function 

The estimation of higher even moments of particle accelerations requires specific data 
processing, as we show in the following. We first discuss the raw distributions of the 
acceleration. In figure IH they are shown for Rex = 815; the particles accelerations have 
been normalized to their variance (whose behavior has been discussed above). 
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To leading order, the distribution functions are v ery similar, as observed in the wind 
tunnel measurements by iQureshi et ai\ (|2007l . l2008l ). There is no reduction to Gaussian 
statistics as the particle size grows well into the inertial range (see also iGasteuil ()2009i) 
with particles on integral size). However, as can be observed in the figure, there is a 
significant infiuence of the flow generation process on the form of the probability dis- 
tribution. In figure 01 the PDF for the smallest particles is ide ntical to that obs erved 
for tracers in the VK flo w opera t ed by Bodenschatz's g r oup, c. f. IVoth et al. and 
in numerical simul ations YeunS (20021); Mordant et al. ( 200461). It is different fr o m the 



PDFs reported by iQureshi et al\ (120071 . 120081 ) and by lAvvalasomavaiula et all (|2006l) 



from wind-tunnel measurements. These observed variation are more pronounced than 
what could be expected from Reynolds number variations alone between the experiments, 
and we thus propose that the PDFs of accelerations are n ot universal, but flow-d ependent . 
This conclusion is linked to the observation (reported in Ouellette et al. ( 20061 )) that the 
Lagrangian small scale dynamics still reflects the anisotropy of the larges scales. 

Investigating the possibility of sub-leading changes in the statistics of acceleration with 
size or Reynolds number can be done by studying higher order moments, starting with 
the distribution flatne ss. It requires a converged measurement of the PDFs and, as shown 
for tracer particles bv lMordant et al\ ([2004a,), this implies extremely large data sets. As 
a flrst attempt, we flt the acceleration PDF with a model functional form which we then 
use to estimate the flatness of the distribution. The procedure is as follows: we assume 
that the statistics is described by a functional form Fe{a); {9} is a set of adjustable 
parameters which are determined by minimizing the distance a^n(a) — a'^J-e{a), where 
n(a) is the measured distribution. Two trial distributions have been tested: 




J-^i^(a) = i_crf 
4^/3 \ 

which stems from the assumption that the acceleration amplitude has a lognormal 
distribution, and the stretched exponential functional form 
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Figure 5. (left): lognormal and stretch exponential fits 11(2; — a/arms) and its moments x^Hix) 
and a;*n(a;) , for data at Re\ = 815 processed using a gaussian smoothing with width w — 20. (o) 

filtered experimental data, ( — ) stretch exponential fit (with flatness F = 55), ( ) lognormal 

fit with (s = 0.96 and F — 68). (right, top): evolution of the flatness estimated from the fitting 
functions, as a function of w, and corresponding linear interpolations to zero width, (right, 
bottom): relative evolution of the two estimators. 



^r(a) = Aexp( , °' I 

(A being a normalization constant) which has three adjustable parameters (a,/3,7), 
and allows for a finer adjustment of the distribution in the tails. Note that with distri- 
butions having such extended wings, a 'brute force' measurement of the flatness factor 
within a 5% accuracy would mean a resolution of the distribution up to about 100 
standard deviations, and events with probability below 10~^^ - clearly outside of direct 
experimental reach! 

Figure [5] shows a comparison of the acceleration pdf for the smallest particles at Re\ = 
815 and the corresponding fits that minimize the distance to the quantity a^n(a). As one 
can see, both functional form fit correctly the experimental data up to aj arms ^ 20, the 
stretch exponential form showing a better agreement with the second order moment. As 
reported in previous studies, the moments of the acceleration PDFs strongly depend on 
the width w of the smoothing kernel used to extract the velocity data from the modulated 
optical signal (cf. section 2). We thus estimate the flatness by fitting the different PDFs 
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Figure 6. Variation of the estimated flatness as a function of D/jy for different sizes and 
different reynolds numbers, (a) stretch exponential estimator; (b) lognormal estimator. 



for decreasing w and then interpolate to zero width. The result of this procedure is 
shown in figure Eljd) for the two trial functional forms. As can be seen the flatness 
derived from the lognormal estimator is roughly 1.5 higher than the one computed from 
the stretch exponential estimator. Indeed, we find that the flatness computed from each 
function follow an identical behavior: when plotting, as in figure[5][e), the values estimated 
from the lognormal fit versus the ones computed from the stretch exponential estimator, 
one observe a linear dependance. The values of the flatness reported are understood as 
estimates, the true values strongly depending on the real shape of the acceleration PDF 
in the far tails. We propose however that the variations detected here yield a first order 
approximation of the evolution of acceleration statistics with particle size. 

The results are shown in Figure |6lja,b) . For both estimators, one observes a reduction 
of the flatness with increasing particle size. The observed decrease can be understood if 
one takes into account the intermittency of the pressure increments for in er tial range sep 
arations. Indeed, following the approach developed bv lVoth et all (|2002l ): lOureshi et al 



(|2007t ). one estimates the acceleration flatness F(D) by assuming that the force acting on 
the particles is dominated by the pressure gradient at size D. All moments of the acceler- 
ation ((a^)) should then scale as Sp{D)/D", its behavior being dictated by the pressure 
structure functions. Now, in order to estimate the pressure structure functions, one can 
either use the ansatz that pressure increments scale as the square of velocity increments, 
(SdP) {{Sdv)'^) or one can measure directly the scaling of pressure in the experiment. 
In the first case, one obtains Sp{D) cx _D^^", C„ being the structure function exponents 
of the velocity increments. One then has F{D) cx Z)'''*"^''* ~ £)-o.42 ^^^^ assumes 



a logn ormal scaling for the Eulerian velocity structure functions as in iChevillard et 



independent of the Reynolds number. The direct experimental measurement of 
pressure (using piezoelectric sensor mounted flush with the lateral wall, in the mid-plane 
of the flow) yields Sj, oc i:»i-2±o.i and oc D-o.38±o.03^ va lues which are in 

agreem ent with Eulerian DNS data at Re\ = 180 (from the DNS by ICalzavarini et al 



] n the case of the stretch exponential estimator, these predictions for the scaling 
exponent are consistent with the the value a ^ —0.4 obtained by fitting the data with a 
power law Fs{D) = A{D/ri)" with D/t] in the range [10,40]. In the case of the lognor- 
mal estimator one finds a scaling law Fi{D) cx {D /r])~^-^ . Here we stress that F}^^{a) 
and J-f^{a) are intrinsically different distributions, therefore it is in principle impossi- 
ble to fit both curves with the same scaling exponent. The true value of the exponent 
(if it exists) should depend on the real shape of the PDFs. However, the consistency 
between the estimation using the stretch exponential estimator and the Eulerian mea- 
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surements of pressure suggests that this estimator is a good fit of the acceleration PDFs 
in the case of large acceleration flatness (F > 20). This was confirmed by comparing 
the quality of the two d ifferent estimators with data obtained at Re\ ~ 180 {F ~ 27.5, 
Calzavarini et al. ( 20091 )). Using a troncated datasct as a test, the stretch exponential 



estimator showed to be a better fit than the lognormal estimator, and was able to give 
an estimate of the Flatness 15% lower than the converged value computed with the all 
dataset. Finally, we note that following the same approach, one can use the second or- 
der pressure structure function to get a new estimate of the decrease of the acceleration 
variance (a|,). Assuming a lognormal scaling for the velocity increments, one then finds 
{a\)) (X ({5dP/D)^) ^ D'^-*"^ = £)-o.78±o.02 ygj,y close to the experimental measure- 
ments which yields (aj^) ~ £)-o.8±o.i^ These two values are in a very good agreement 
with the best fit shown as a dashed line in figure 2 which yields (ajj) /{a^) cx {D /rj)^'^-^^ . 
In our opinion this is an indirect proof that intcrmittency should play a role, and that 
our results concerning the acceleration flatness, if preliminary, are consistent. 



4. Concluding remarks 

Scaling properties of the acceleration of (neutrally buoyant) inertial particles, i.e. con- 
cerning finite size effects for the motion of advected particles, have been investigated in 
detail. One finding is statistics is well accounted for by the behavior of pressure incre- 
ments, which in turn arc connected to the Eulerian properties of velocity increments. A 
first estimation concerning the evolution of the particle acceleration PDF with size sug- 
gests a reduction of the fiatness with a (D/rj)"'^-'^ scaling behavior. However, since the 
precise functional form actually depends on the large scale fiow properties, this results 
need to be confronted to measurements in other fiows, such as grid or jet turbulence. 
Another observation is that the evolution of the acceleration auto-correlation functions 
allow for the estimation of a particle response time. Our finding is that it increases lin- 
early with the particle size in the range of sizes explored. One point that will deserve 
further studies is the dependence on Reynolds number. Here we have achieved a satis- 
factory collapse of measured quantities when the particle size is made non-dimensional 
using the Kolmogorov length This may be justified for the R\ values and particle sizes 
investigates here; for larger particles or lower Reynolds numbers one may also want to 
probe scaling properties as a function oi D/L, with L the flow integral scale. Here, L was 
much larger than the largest particles size, but such may not be the case for all studies. 
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